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We analytically demonstrate helicity determined selection rules for intershell tunneling in double- 
walled nanotubes with commensurate (c-DWNTs) and incommensurate (i-DWNTs) shells. For 
i-DWNTs the coupling is negligible between lowest energy subbands, but it becomes important as 
the higher subbands become populated. In turn the elastic mean free path of i-DWNTs is reduced 
for increasing energy, with additional suppression at subband onsets. At low energies, a Luttinger 
liquid theory for DWNTs with metallic shells is derived. Interaction effects are more pronounced in 
i-DWNTs. 
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Due to their unusual physical properties, cf. e.g. Q 
carbon nanotubes have attracted lots of attention. Car- 
bon nanotubes can be single-walled (SWNT) or multi- 
walled (MWNT), depending on whether they consist of a 
single or of several graphene sheets wrapped onto coaxial 
cylinders, respectively. Electronic properties of SWNTs 
are mostly understood [lj. In particular, SWNTs are 
usually ballistic conductors 0], and whether a SWNT is 
metallic or semiconducting is solely determined by its so 
called chiral indices (n, m). Due to the one-dimensional 
character of the electronic bands at low energies, Lut- 
tinger liquid features at low energies have been predicted 
0,0] and observed 0, 0. The situation however, is 
much less clear for MWNTs. Except for few experi- 
ments, see e.g. 0, MWNTs are typically diffusive 
conductors, see e.g. 

ana, 

with current being carried 
by the outermost shell at low bias [ToL Hl| and also by 
inner shells at high bias E3 • Intershell conductance mea- 
surements consistent with tunneling through orbitals of 
nearby shells have recently been reported EH- More- 
over, which kind of electron-electron correlation effects 
determine the observed zero-bias anomalies [3 Hi of 
MWNTs is still under debate. To better understand 
these features, i.e., the role of mierehell coupling on 
transport properties of MWNTs, some experimental E3 
and theoretical El H H |H HJ H works focussed 
on the simplest MWNT's realization, namely on double- 
walled nanotubes (DWNTs). One main outcome is that a 
relation must exist between the intershell coupling, shell 
helicity and transport properties. Specifically, two shell 
are called commensurate (incommensurate), if the ratio 
between their respective unit cell lengths along the tube 
axis, is rational (irrational) UjL For example, using tight- 
binding models, Saito et al. |17| numerically found energy 
gaps opened by the intershell coupling in a DWNT with 
two armchair (and hence commensurate) shells. Ab-initio 
calculations confirmed these results. In general, 

numerical evidence of a negligible intershell coupling in 
DWNTs with incommensurate shells (i-DWNT) at low 
energies is found 

[ni ei mini. 



In this Letter, we derive an analytical expression, yield- 
ing helicity- dependent selection rules for tunneling, for 
the effective intershell coupling. For i-DWNTs the inter- 
shell coupling is negligible between the lowest subbands 
but it becomes important when higher subbands are in- 
volved. We show that this in turn yields an elastic mean 
free-path which decreases with energy and which shows 
a characteristic suppression at subbands onset. Then by 
including intra- and inter-shell Coulomb interactions, we 
show that metallic DWNTs can be described by Lut- 
tinger liquid theory at low energies. The tunneling den- 
sity of states has a power-law behavior with different ex- 
ponents for i-DWNTs and commensurate-shells DWNTs 
(c-DWNTs). 
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FIG. 1: (a) Graphene lattice. The x and y axes are along the 
armchair and zigzag axes respectively. The distance between 
two nearest carbon atoms is a<j ~ 1.42 A. The unit lattice 
vectors are ai and a2. The vectors di connect three nearest 
neighbour atoms, while p and r are two vectors required to 
specify the position of a carbon atom, (b) Cross section of 
a DWNT. Atoms A and B in two shells of radii R a and Rb, 
respectively, are projected onto this cross section. 

To derive the helicity-dependent selection rules, we use 
a tight-binding model for non interacting electrons with 
one 7r-orbital per carbon atom EI and follow |24|. This 
model is described by the Hamiltonian 



(ij) «' 



bj C ai C bj 



H.c. 



(1) 



where (3 = a, b is the shell index, (ij) is a sum over 
nearest neighbors in a shell, 70 ~ 2.7 eV 1] is the in- 
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trashell nearest neighbour coupling. The intershell cou- 
pling is U ai<rbj = t e~ [d( - r <"- r i. 3 )- A )/ Q * 7 where t ~ 1.1 eV, 
A ~ 0.34 nm, d(r a i , ?b±L is the distance between two 
atoms, and at ~ 0.5 A [l|. We introduce the transfor- 
mation 

V k 



tS(q va R a — q v bRb)S{q ua — (fob), 



where r\ = ± is the index for the two interpenetrating 
sublattices in a graphene sheet, and Np is the number of 
graphene unit cells on a shell. The Hamiltonian reads 



Ha = X^ k 4»k c /3-^k 

/3r,k 



J f)i?bk 6 



7o e zkd J = 



where the mirashell coupling is 7k = X) 
|7k|s, with dj the vectors connecting the three nearest 
neighbour carbon atoms. Introducing vp = R + X^, 
with R a lattice vector, = Pp + fjpT, where p a — p b 
describes the relative position of the two shells, cf. Fig. 
1, the elements of the mierehell 2x2 coupling matrix are 

T VaVb (k a ,k b ) = £ e* G - x -* G - x ^ ka+GQ , kb+Gb . (4) 

G a G(, 

Here G is the graphene reciprocal lattice vector 
4tt 



G 



3ao V 2 

with h,l 2 = 0,±1,±2, .... Finally, 



(h-k), \(h+h)) , 
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with A ce n the area of a graphene unit cell. We no- 
tice that the Hamiltonian Eq. (|2J) is diagonalized, in 
the absence of intershell coupling, by the transformation 
U = /* ). For later purposes, cf. Eq. Q, we call 

^v a v b = (U*TU) Valli the elements of the intershell tun- 
neling matrix between two Bloch states in different shells. 
Here v — =p is the index for bonding/ antibonding states 
corresponding to negative/positive energies ep^(k) with 
(3 = a, b, respectively. 

It is convenient to introduce coordinates u and v, which 
are along the tube axis and the circumference direction 
respectively, cf. Fig^b). Then k a • v = k va obeys £ a — 
kvaRa, due to the periodic boundary conditions along the 
circumference. Likewise, £ b — k vb R b . Here, the integers 
£ a and lb characterize energy subbands. In contrast, k u = 
k • u is continuous, cf. Fig.|2£b). The distance between 
two atoms A — {R a cos(v a / R a ), R a sm(v a /R a ), u a ) and 
B = (R b cos(v b /R b ),R b sm(v b /R b ),u b ) is then d(r a , r b ) = 
D(v a /R a - v b /R b ,u a - u b ) where 



D(zi,z 2 ) 



R a -R b \ 2 + AR a R b sin 2 ( Zl /2) + {z 2 f 



We then find iq a , q! 
where the prefactor t is calculated as 

'dz x dz 2 e -(D(z^)-&)/a t 



t = tn 



cell 



Therefore, according to the two S- functions in iq a , qi> , the 
effective intershell couplings 7^ a% (k a ,kt) between two 
shells (n a ,m a ) and (n b ,m b ) are nonzero if they satisfy 
the following selection rules, 

L + (nJia + m a l 2a ) = £ b + {ribhb + m b l 2b ), (5a) 
k ua + T{n a ,m a ) = k ub + T{n b ,m b ), (5b) 



H.c. , (3) with T(n, m) 



3aoC(n,m) 

v3ao£(n,m), with C(n,m) = y/n 2 + 



((n+2m)ii — (2n+m)l 2 \. Here 



nm, is the 



circumferential length of shell (n,m). At low energies 
only the lowest subband determined by Z£p = 2np + mp 
in each shell is important, which fixes the values £ a and £ b - 
For c-DWNTs, e.g. if the two shells are either both arm- 
chair or zig-zag, Eqs. (JSJ can always be satisfied. More- 
over, the dominant contribution is for k ua = k ub . On 
the other hand, for i-DWNTs, e.g. a (9,0)@(10, 10), the 
selection rule Eq. (|5b|) can only be satisfied if the dif- 
ference k ua — k ub takes finite values. At low energies, 
this condition is never met. At higher energies higher 
subbands must be considered as well, and the selection 
rules can be satisfied. Notice that whenever the l.h.s. 
and r.h.s. of Eqs. JSJ are not close to zero, the effec- 
tive intershell coupling is exponentially suppressed |25| . 
We show now that, for i-DWNTs, the increase of the 
inter-shell tunneling is at the origin of an elastic mean- 
free path l e i which decreases with increasing energy, and 
which shows a characteristic suppression at each sub- 
band onset. Our analytical results are in agreement with 
recent ab-initio calculations, showing a cross-over from 
ballistic to diffusive behavior in DWNT as the energy 
increases [13, EI , as well as with the experimental ob- 
servation that MWNT mostly exhibit diffusive behavior. 
To evaluate the elastic mean-free path l e i,b(E) = VFTb{E) 
for electrons in the shell b, the life-time r b (E) for elec- 
trons with energy E is needed. Here vp = 8 x 10 5 m/s 
is the Fermi velocity for nanotubes. To be definite, 
h/ n (E) = J2^=±W T bM,)S{E - £ M (k)), with e 6i „(k) 
the dispersion relation in shell b, and 

-= E [ dk' u \T^,^,(e b Ak))\ 25 (ebA^-^(k')). 



Tfe,k 



' = ± 



(6) 

For DWNTs the 4 x 4 T-matrix [2j| is evaluated to be 

Tkk'M = v(k,k') + E v ( k ' k i) G ( k i^) T k lk 'M, 



where V^ b (k,k') = (V^(k,k / ))* = % aVh {^\d) 
and V^, (k, k') = 0. The elements of the retarded 
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Green's function G are G^,(u>, k) = (w — ep iU (k) + 
iO + )~ 1 5p l 3i5v U i . For i-DWNTs T VaUb couples subbands 
with different energies, cf. Fig. 2. Thus, in general is 
£6,i, (k) ^ e a , V '{k'), i.e., (5(e b ^(k) - £ QjI/ /(k / )) = 0, and 
f3 = b in Eq. (JSJ). Hence, to lowest order in T, and if 
£&,„(k) 7^ £ 0) ^/(k'), the life-time Tb^ v is obtained insert- 
ing in Eq. © 



. ^ , £6,r/(k) -£„,„! (ki) + ^ 



(7) 

The elastic mean free paths i e i, (10,10) ancl ^ei,(9,o) f° r 
electrons in the outer and inner shell, respectively, of a 
(9,0)0(10,10) DWNT are shown in Fig. 3. It is clearly 
shown that before the first subband onset, the motion is 
ballistic also for i-DWNTs of lengths up to ~ 5pm. 

In the remaining of the paper we consider only metal- 
lic shells and include electron-electron correlation effects 
in the low energy regime where only the first subband 
of each shell is populated. In this regime transport is 
ballistic for c-DWNTs as well as for i-DWNTs. Due to 
the linearity of the dispersion relation, a multi-channel 
Luttinger liquid description can be used |2(| [2^ ■ 
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FIG. 2: (a) Schematics of the effect of the selection rules for 
i-DWNT (upper) and c-DWNT (down), (b) Energy subbands 
of the i-DWNT (9, 0)@(10, 10) in the absence of intershell cou- 
pling. The effect of the latter is to induce transitions between 
different subbands in different shells. The dominant coupling 
of the lowest armchair (zig-zag) subband is indicated. 

At first, we consider an i-DWNT where the intershell 
coupling can be ignored. The unperturbed Hamiltonian 
can be written as 



H = -ihv F du ^L*f3 d u*Pr a *t3 



(8) 



where r = ± is the index for right/left movers, a = ± for 
the two independent Fermi points of a shell, and a = ± 
for up and down spins. The electron density operator 
is pulu) = E™<T/3^a CT /3( u )V' ra ^(«)- The two shells are 
only coupled by the Coulomb interaction, which gives 
rise to forward, backward, and Umklapp scattering pro- 
cesses. Experimentally, the Fermi points of nanotubes 
are usually shifted away from the half-filling due to dop- 
ing or external gates. We assume that this is the case and 
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FIG. 3: Elastic mean free paths for electrons in the outer and 
inner (inset (a)) shells of the i-DWNT (9, 0)@(10, 10). Notice 
the dips in correspondence of the first two subband onsets. 



hence neglect Umklapp processes. Since we are not inter- 
ested in the extremely low temperature case, the back- 
ward scattering processes are also ignored here [2^. In 
the following, we only consider forward scattering pro- 
cesses described by the Hamiltonian 



\ ^ I ' dudu' pp(u)V0i3>(u- u')p/3>(u'), 



(9) 



with the effective one dimensional interaction 
f 2 ^ f 2 ^' dvdv' 

where U(r) is the Coulomb interaction. The Hamiltonian 
H = H + Hps can be diagonalized by the bosoniza- 
tion procedure discussed in Ref. [j^. We introduce 
bosonic field operators for the total/relative (S = ±) 
charge/spin (j = c,s) modes in shell f3, as well the to- 
tal/relative (£ — ±) modes with respect to the two shells 
obeying the commutation relation [9j^(u), 4>j>8'£'(u')] — 
— (i/2)Sjj'Sss'S^isgii(u — u / ). The Hamiltonian H can be 
then decoupled into 8 modes as 



J2^r i du ( K m{ d v®m( u )) 2 +j^{du<t>m{ u )) 2 

M J \ 

(10) 

Only the two total charge modes are renormalized by the 
Coulomb interactions with velocities v c +± = vf/Kc+± 
and interaction parameters 



( n ) 

where V/3/3' = V/3{3'(2ir/L) is the Fourier transform of the 
interaction potential at long- wave lengths, with L the 
nanotube length. The remaining modes are neutral with 
parameters VjS£ — vp and Kj$£ = 1. 
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We consider now c-DWNTs where intershell tunneling 
is relevant. The intershell tunneling Hamiltonian is 



Ht = 



J2^o J du ipl aaa (u) ip raab (u) + H.c 



(12) 



where for simplicity the tunneling element 7++(k, k') 
is evaluated at k = k' = K with the Fermi point K 
of graphene, and is the constant %. As detailed in 
|2fjj . the Hamiltonian Ha + H t can be exactly diago- 
nalized. One finds the same form as in Eq. JHJ) where 
now the index (3 — 0, tt stands for bonding and anti- 
bonding states, respectively. Moreover, the Fermi wave 
vectors of the two independent Fermi points are shifted 
as kp — ► k F ± (Tq/Hvp) where ± stand for tt and 0, 
respectively. We retain again only (intraband and in- 
terband) forward scattering described by the Hamilto- 
nian Eq. JjJ}, where now the scattering potentials Vppi 
are V m = VW = V „/2 = (V aa + V bb + V a b)/i, so that 
bosonization brings again the total Hamiltonian in the 
form Eq. with 6 neutral modes and 2 renormal- 

ized total charge modes. The tunneling density of states 
(TDOS) of both shells, Pb/ a { e ), immediately follows |2flj . 
For i-DWNT is 

Pb/ a (e) ~ \e\ a »'° , 

with exponents a b / a being different for electrons tun- 
nelling into the middle or end of a nanotube: 



ad,6/a- lE^^" 1 



Oibu\kM/a = g ^ A b / a ^K c+ £ + — - 2 



(13) 



Here the coefficients A„ , A b = 1 — A a are related to the 
eigenvalue problem [2al27| . For c-DWNT is 

p b/a {s)~\e\ a ° + | e |^e(e-2r ), 

where O(x) is the Heaviside step function and 27o is the 
gap between antibonding and bonding states. Because 
the intraband forward scattering potentials are equal, 



is ao = a-K- We find a e 



id/bulk 



given by Eq. 1 (15)1 with 



A b = A a = 1/2. For illustration we calculate the tunnel- 
ing exponents for the (10,10) shell of a (9, 0)@(10, 10) 
and of a (5,5)0(10,10) with radii R a « 3.4 A and 
R b ~ 6.8 A. We find a on d = 1.21,abuik = 0.50 for a 
(9, 0)@(10, 10) DWNT and a cnd = 0.80, a bu ik = 0.34 for 
a (5,5)@(10, 10) DWNT. For comparison, for a (10,10) 
SWNT is a on d = 1-25 and abuik = 0.52. Hence, the ex- 
ponents of DWNTs decrease due to the screening effect 
of the inner shell with respect to a SWNT. The intershell 
coupling reduces the exponents further. Notice that for 
Fermi liquids is a cn d/buik = 0. 

In summary, we derived selection rules according to 
which the intershell coupling is only negligible in i- 
DWNTs at low energies. An analytical expression in 



Born-approximation for the elastic mean free path was 
provided. Including the Coulomb interaction, we de- 
veloped a low energy Luttinger liquid theory for metal- 
lic DWNTs according to which the intershell coupling 
strongly reduces the tunneling density of state exponents 
in c-DWNT with respect to those of i-DWNTs. 

The authors would like to thank G. Cuniberti, J. Keller 
and C. Strunk for useful discussions. 
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